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ABSTRACT 

The AKARI Deep Field South (ADF-S) is a ~12 sq. deg. region near the South 
Ecliptic Pole that has been observed with deep scans in the far-infrared by the AKARI 
satellite. As such it is becoming one of the key extragalactic survey fields. We here 
present complementary observations of the ADF-S conducted by the Spitzer Space 
Telescope at wavelengths of 24 and 70 /im. We extract source catalogs at each of 
these wavelengths reaching depths of — 0.2mJy at 24/im and — 20mJy at 70/im. We 
also apply an K-to-24/im colour criterion to select objects with galaxy-like colours in 
the 24/im survey. Completeness corrections as a function of flux density are derived 
for both catalogs by injecting artificial sources of known flux density into the maps, 
and we find that our surbveys are 50% complete at 0.26mJy and 24mJy at 24 and 
70/im respectively. We can thus produce number counts as a function of flux density 
for the ADF-S at 24 and 70/im. These are combined with existing literature counts 
and compared to four different number count models derived from galaxy evolution 
models. One complicating factor for the ADF-S counts is the presence of a foreground 
galaxy cluster at z=0.04 in the field. We examine the ranges of flux densities to which 
this cluster might make a contribution to the counts and find hints that the 24/im 
luminosity function of the cluster galaxies might be enhanced above that of field 
galaxies. Full catalogs for these ADF-S Spitzer surveys at 24 and 70/im are made 
available as part of this paper. 
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1 INTRODUCTION 

The Japanese AKARI mission (Murakami et al., 2007) has 
carried out large legacy observations in high visibility ar- 
eas at the North and South Ecliptic Poles. At the South 
Ecliptic Pole (SEP), AKARI has made a deep survey at 
far-infrared wavelengths with its Far-Infrared Surveyor in- 
strument (FIS, Kawada et al., 2007) at wavelengths of 65, 
90, 140 & 160/tm to 5a sensitivities of 46, 15, 180 & 600mJy 
respectively (Shirahata et al., 2009, PI: S. Matsuura). The 
AKARI SEP survey, hereafter referred to as the AKARI 
Deep Field South (ADF-S) covers approximately 12 square 
degrees centred at R.A.=4 h 44 m 00 s , Dec=-53deg 20'00."0 
(J2000) (Matsuura et al., 2010), a region where the cir- 
rus brightness is spectacularly low < Bioo >~ 0.2 MJy/sr 
(Schlegel et al., 1998). The main science drivers for the sur- 



vey are the analysis of clustering and large scale structure 
using the large contiguous area, and fluctuation analysis of 
the background. In addition to the far-infrared data, shal- 
low mid-infrared data also exists over the entire area, along 
with deep optical _R-band imaging to 25th magnitude. The 
entire ADF-S has also been covered by the Balloon-borne 
Large- Aperture Submillimeter Telescope (BLAST, Pascale 
et al. (2008) at 250, 350 & 500/mi. Other ancillary data 
on smaller scales include a central area of approximately 
1 square degree covered by the AKARI Infra-Red Camera 
(IRC, Onaka et al., 2007) at wavelengths of 3, 4, 7, 11, 
15 &24 /mi, optical UBVI imaging, submillimetre observa- 
tions with the Atacama Pathfinder Experiment (APEX) us- 
ing the Large Bolometer Camera (LABOCA), the Atacama 
Submilimetre Telescope Experiment (ASTE) using the As- 
tronomical Thermal Emission Camera (AzTEC), and 20cm 
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Australia Telescope Compact Array (ATCA) radio observa- 
tions (White et al., 2009). 

Spitzer Space Telescope (Werner et al., 2004) obser- 
vations of this field, covering almost the entire 12 square 
degrees, were also conducted. These observations were car- 
ried out with the Multiband Imaging Photometer for Spitzer 
(MIPS, Rieke et al., 2004) instrument at wavelengths of 24 
and 70/im. Observations in the 160^im MIPS channel were 
not collected since the AKARI data at 140/im is at compa- 
rable depth. We here present our processing of the Spitzer 
ADF-S data and the catalogs extracted from our reduced 
images. The catalogs are made available for use by the astro- 
nomical community to maximize the scientific return from 
this Spitzer data and from the ADF-S in general. 

As well as producing the Spitzer images and catalogs 
for the ADF-S, we also extract number counts at 24 and 
70/mi for this field and use them to set extra constraints 
on models of the galaxy population at these wavelengths. 
Previous work on counts at 24/im includes both deep counts 
from small areas (Papovich et al., 2004; Le Floc'h et al. 
2009) and shallower counts from the largest Spitzer sur- 
vey, the Spitzer Wide-Area Infrared Extragalactic Survey 
(SWIRE, Shupe et al., 2008) while Bethermin et al. (2010) 
present a compendium of counts at 24/im over a wide range 
of fluxes, alongside 70 and 160 /jm counts. The observations 
we present here are reasonably well matched to the SWIRE 
survey counts in terms of depth, going somewhat deeper 
but covering only one quarter of the SWIRE area. Previous 
70/mi counts have concentrated more on small area fairly 
deep surveys, including COSMOS (Frayer et al., 2009), and 
GOODS-N (Frayer et al., 2006a). The largest area survey 
so far fully published is the FLS field (Frayer et al., 2006b) 
which covers 4 square degrees to depths comparable to our 
ADF-S Spitzer data. The SWIRE survey counts at 70/im 
have yet to be fully analyzed and published, though prelim- 
inary results for high completeness (S/N >10) SWIRE 70/im 
sources are presented by Frayer et al. (2009). Our work thus 
provides useful additional constraints for both the 24 and 
70/im counts largely at bright and intermediate flux densi- 
ties. 

One problem that arises in the ADF-S field which has 
an impact on counts studies is the presence of a foreground 
rich galaxy cluster in the field, DC 0428-53 (also known as 
Abell S0463) (Dressier, 1980; Dressier & Schectman, 1988) 
at a redshift of ~ 0.04. Galaxies in this cluster will tend 
to enhance the brighter counts, and we observe this effect. 
Any large area survey will be subject to the sample variance 
induced by large scale structures, so this is not a unique 
problem. Nevertheless, the presence of this cluster affects 
our interpretation of the brightest counts. It also yields the 
opportunity to study known cluster members in the infrared 
- something we plan to do in future publications. 

The remainder of this paper is organized as follows. In 
Section 2 we describe the Spitzer observations and the data 
processing steps used to produce the reduced maps. In Sec- 
tion 3 we describe source extraction at 24/im and the pro- 
duction of the 24/im counts. Section 4 does the same at 
70/im. Section 5 describes the source catalogs produced us- 
ing these data, while our results are discussed and the counts 
compared to a variety of models in Section 6. In Section 7 we 
draw our conclusions. Throughout this paper we assume the 



concordance cosmology with Ho = 72kms Mpc 1 , Q. m = 
0.3 and ft A = 0.7. 



2 SPITZER OBSERVATIONS AND DATA 
PROCESSING 

The Spitzer observations were performed using the MIPS 
scan map mode under programme ID 50581 PI: Mark De- 
vlin and consisted of a total of 91. 2h of observing time. 
The instrument can simultaneously acquire data in all three 
MIPS bands (24, 70, and 160 /im). Since only the 24 and 
70 /im data were needed for comparison to the AKARI 
data, the instrument was used in a cryogen-conserving mode 
whereby the telescope warmed up slightly to a point where 
the 160 /im data are not usable. The observations consist of 
a total of 34 Astronomical Observation Requests (AORs). 32 
of the AORs were scan map observations with 9 scan legs 
that were 1.5° long and that were offset from each other 
by 160 arcsec. The other 2 AORs were scan map observa- 
tions with 5 scan legs that were 0.5° long and that were 
offset from each other by 160 arcsec. All scan maps were 
performed at the medium scan rate, and all scans were ap- 
proximately parallel to the short axis of the SEP field. The 
observed field is approximately 6.5° x 2.2° and contains 11.8 
square degrees of data at 24 /im and 11.5 square degrees of 
data at 70 /im. Execution of these AORs was spread over the 
period from 2008 September 24 to 2008 October 1. Typical 
total exposure time per map pixel after these observations 
were completed was 164s at 24/im and 73.7s at 70/im. 

The 24 and 70 micron images were created from raw 
data frames using the MIPS Data Analysis Tools (DAT) 
version 3.10 (Gordon et al. 2005) along with additional pro- 
cessing steps. The data processing for the two wave bands 
differs significantly, so the processing for each band is de- 
scribed in separate paragraphs below. 

The individual 24 /im frames were first processed 
through a droop correction (removing an excess signal in 
each pixel that is proportional to the signal in the entire 
array) and were corrected for non-linearity in the ramps. 
The dark current was then subtracted. Next, scan-mirror- 
position dependent flats were created from the data in each 
AOR and were applied to the data. Detector pixels that 
had measured signals of 2500 DN s _1 in any frame were 
masked out in the following three frames so as to avoid hav- 
ing latent images appear in the data. Next, a scan-mirror- 
position independent flat was created from the data in each 
AOR and were applied to the data. Following this, we sub- 
tracted variations in the background related to the position 
of individual frames in the scan legs, which removes zodiacal 
light emission and variations in the background related to 
the scan mirror position as well as any other cyclical back- 
ground fluctuations. Next, a robust statistical analysis was 
applied in which the values of cospatial pixels from different 
frames were compared to each other and statistical outliers 
(e.g. probable cosmic rays) were masked out. After this, a 
final mosaic was made with pixel sizes of 4 arcsec 2 , and the 
data were multiplied by the flux density calibration factor 
(4.54 ±0.18) x 10~ 2 MJy sr" 1 [MIPS instrumental unit]" 1 
(Engelbracht et al., 2007). 

At the beginning of the 70 /im data processing, ramps 
were fit to the reads to derive slopes. In this step, readout 
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jumps and cosmic ray hits were also removed, and an elec- 
tronic nonlinearity correction was applied. Next, the stim 
flash frames (frames of data in which a calibration light 
source was flashed at the detectors) were used as responsiv- 
ity corrections. After this, the dark current was subtracted 
from the data, and an illumination correction was applied. 
Following this, short term variations in the signal (often re- 
ferred to as drift) were removed, and additional periodic 
variations in the background related to the stim flash cy- 
cle were subtracted; this also subtracted the background 
from the data. Next, a robust statistical analysis was applied 
to cospatial pixels from different frames in which statistical 
outliers (which could be pixels affected by cosmic rays) were 
masked out. Once this was done, final mosaics were made 
using pixel sizes of 4 arcsec 2 . Finally, the data were mul- 
tiplied by the flux density calibration factor 702 ± 35 MJy 
sr~l [MIPS instrumental unit] -1 (Gordon ct al., 2007). 

The final 24/im image together with a coverage map is 
shown in Fig. I, while Fig. 2 shows the final 70/im image 
and coverage map. A magnified section of the image at each 
wavelength is shown in Fig. 3 where it is easier to see the 
full range of fluxes of the objects detected in the survey. 



3 24 MICRON SOURCE EXTRACTION AND 
COUNTS 

3.1 Source Detection and Extraction 

For the detection and extraction of sources at 24/im, we fol- 
low the recipe given in Shupc ct al. (2008) since our data 
is broadly similar to the 24/tm SWIRE survey data dealt 
with in that work. Source extraction was performed using 
SExtractor (Bertin & Arnouts, 1996). We use a local back- 
ground calculated over a region of size 128 x 128 pixels. 
Sources were detected using a 2a threshold with a minimum 
five connected pixels above this threshold. Aperture pho- 
tometry using a 5.25 arcsec aperture and a Kron flux (Kron 
1980) were extracted for each source. The flux assigned to a 
source was the aperture flux unless the measured size of the 
source was >100 pixels and the SExtractor stellarity prob- 
ability was less than 80%. If these were both true the Kron 
flux was used instead of the aperture flux, as in Shupe et 
al. (2008). We only extract sources from regions covered by 
at least ten MIPS frames. This results in the detection of 
41503 sources in a region of 11.26 sq. deg. 

We follow the SWIRE team in applying a correction fac- 
tor of 1.15 to match the MIPS team calibrations to bright 
stars (Surace et al., 2005). We also adopt their aperture 
correction factor of 1.55 to account for light outside the 
point source 5.25 arcsec aperture, a factor of 0.961 to correct 
the photometry to an assumed F v ~ source spectrum 
(Stansberry et al., 2007) and an additional correction factor 
of 1.018 to account for changes in the 24/im MIPS calibra- 
tion factor (Surace et al., 2005; Shupe et al., 2008; MIPS 
Data Handbook, 2006)). 

3.2 Star-Galaxy Separation 

At 24/im flux densities brighter than a few mjy stars are 
a potentially serious contaminant for investigations of ex- 
tragalactic number counts (Shupe et al., 2008). To separate 
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Figure 4. Colour-Magnitude Diagram of 2MASS K band and [24] 
for the ADF-S field. The dashed line shows the colour cut used to 
separate out the stellar contribution to the number counts (below 
the line). 

stars from galaxies, we cross match the sources in our 24/im 
catalog to sources in the 2MASS K band catalog (Skrutskic 
et al., 2006) using a 3 arcsecond search radius. Of a total of 
41503 24/im sources extracted in areas covered more than 
10 times by MIPS, we obtain K band matches to 2206. The 
resulting K-[24] colour magnitude diagram is shown in Fig- 
ure 4, where we adopt [24] as a 24/im magnitude with a zero 
point of 7.43 Jy in colour corrected units. We use a colour 
cut of K-[24]< 2, as shown in Figure 4 and similar to Shupe 
et al. (2008), to exclude stars. This leaves us with 40556 
sources with galaxy-like colours, and 947 objects with star- 
like colours. The latter are excluded from the main number 
count analysis. Differential counts for the sources identified 
as stars are shown in the Appendix. 



3.3 Completeness Correction and Reliability Test 

The completeness of our 24/im catalog is assessed by adding 
artificial sources of known flux to the map and then deter- 
mining the fraction of these sources recovered by the ex- 
traction process. A thousand artificial sources are added at 
random positions at each of nine different flux density values 
using an empirical noiseless point spread function (PSF) im- 
age. This empirical PSF was created using archival Spitzer 
data of 3C 273, 3C 279, and BL Lac. Data from individual 
AORs for these three galaxies were processed into images 
using the MIPS DAT and then combined so as to filter out 
any extended, asymmetric emission. See Young et al. (2009) 
for additional details. Sources were then extracted from the 
map as described above, and we determined what fraction of 
the input sources were recovered. The resulting completeness 
correction as a function of flux is shown in Fig. 5. We use 
linear interpolation from this data to calculate completeness 
corrections for the differential number counts. We find that 
our catalog is complete at the 80% and 50% levels at fluxes 
of 0.32 and 0.26 mjy respectively. We also use this data to 
check the flux density correction factors we have adopted 
from the SWIRE team ie. the 1.15 correction to bright stars 
and the 1.55 aperture correction factor for a total correc- 
tion factor of 1.78 before colour corrections are included. 
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Figure 1. The ADF-S Field at 24/im. Top: 24^tm image. Bottom: Coverage map at 24/im. Darker colours indicate higher flux density 
or greater coverage respectively. 



We find that this flux density correction factor works well 
at all flux densities except those below ~0.3 mjy where our 
completeness is very poor. 

To assess the reliability of our source extraction we do 
the following. Firstly we subtract the mean per-pixel value 
of the 24pm map from the map. This produces a map with 



mean value zero. We then multiply this map by -1 and add 
back the mean we have subtracted. This produces a map 
with identical noise and statistical properties to the real 
map, but with all the real sources made negative. Sources 
detected by SExtractor in this map are thus false, resulting 
from the noise properties of the map, and the number of such 
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Figure 2. The ADF-S Field at 70/xm. Top: image. Bottom: Coverage map of 70^tm. Darker colours indicate higher flux density or 

greater coverage respectively. 



sources detected can be used to assess the reliability of the 
source catalog. We then process the resulting false source 
catalog in an identical way to the real data, with the ex- 
ception of star galaxy separation since there is no reason to 
expect the false sources to be associated with real objects on 
the sky. We find that our reliability is close to 100% for the 



full flux range used for counts, with our faintest and most 
unreliable flux bin, at 0.285mJy, still having a reliability of 
96%. 
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Figure 3. Enlarged section of ADF-S Field at 24/xm (left) and 70/jm (right) allowing more of the fainter sources to be seen. 
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Figure 5. 24/im Completeness function. Completeness as a func- 
tion of flux density for the 24/rni observations. Completeness is 
determined by injecting artificial sources into the map and then 
determining what fraction are recovered by our source detection 
technique. A thousand artificial sources are used for the determi- 
nation of each point on the completeness function. 



3.4 24/im Counts 

The 24/im differential counts for the ADF-S field are shown 
in Fig. 6, alongside a variety of data (Papovich et al., (2004), 
Shupe et al., (2008)) and models from the literature (La- 
gache et al. (2004), Rowan-Robinson (2009)) as well as the 
Pearson number count models (Pearson & Khan, 2009), de- 
scribed in somewhat more detail below. It should be noted 
that this is not a comprehensive list of available models. 
Many others are available, but we restrict the comparison to 



these models to compare the current 'baseline' model, from 
Lagache et al. (2004) to two more recent models based on 
SWIRE and other recent Spitzer observations while keeping 
the comparison diagrams simple enough for easy interpreta- 
tion. The counts are also given in Table 1. Our data at 24/im 
largely covers flux density regions already covered by other 
surveys, but our greater area results in a larger number of 
sources and thus smaller Poisson error bars. Additional er- 
rors in the determination of the counts can come from large 
scale structures (LSSs) in the galaxy distribution. The clus- 
ter DC0428-43, which lies in our survey region and which is 
discussed more below, is a specific example of this effect, but 
generic LSS will also add to the uncertainty in the counts. 
We assess the uncertainties in the counts coming from LSS 
effects by extracting counts in a number of sub-regions in 
our survey field and performing a counts-in-cells analysis 
(see eg. Wall & Jenkins, 2003). These are then added in 
quadrature to the Poisson errors. The LSS errors dominate 
the Poisson errors in all except the brightest bins. At the 
brightest flux densities, which are not well constrained by 
other surveys, we find a moderate excess of counts beyond 
what is expected on the basis of models. This, however, is a 
result of the foreground cluster at z=0.04 that covers ~l/4 
of this field (see below). In other respects, our 24/xm counts 
are in good agreement with previous observational results. 
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Figure 6. Euclidean Normalized Differential Counts for our 24/tm data compared to other counts in the literature and to theoretical 
models. Data from this work is shown as solid dots for raw counts and open circles for the counts where objects identified with cluster 
galaxies in DC 0428-43 have been removed. Literature counts are shown as follows: open squares: Papovich et al., (2004); open triangles: 
Shupe et al. (2008). The models are plotted as lines, including: solid line: Rowan- Robinson (2009) model; dashed line: Pearson et al. 
(2009) burst mode model; dash-dot line: Pearson et al. (2009) bright mode model; dash-three-dots line: Lagache et al (2004) model. The 
horizontal line above the counts indicates the range of flux densities expected for 0.2L? 2A < L24 < 2-Lj4 objects in the z=0.04 cluster in 
this field and for which some enhancement in counts might be expected (see section 6.3). 



4 70 MICRON SOURCE EXTRACTION AND 
COUNTS 

4.1 Source Detection and Extraction 

Sources are detected in the 70/im map using the APEX im- 
age extraction tool developed by IPAC as part of MOPEX 
for use with MIPS-Ge data (Makovoz & Marleau, 2005) with 
an approach similar to that used by Frayer et al. (2006a, 
2006b, 2009). This uses an empirical point source response 
function (PRF) to extract source flux densities. Sources were 
extracted down to 4cr significance in regions where the sky 
is observed by MIPS at least 10 times. We remove spuri- 
ous low significance sources by cross-matching this 70/im 
source list to the 24/im source list discussed above using a 
9" matching radius. Since the latter is more than 100 times 
deeper in terms of flux density we expect very few genuine 
70/im sources to lack 24/im detections. This process excludes 
404 out of 1702 initial detections, the majority of which are 
in regions of the map with poor 70/im coverage and thus 
enhanced noise. A further 95 sources are excluded because 
there is a small region where the 70/im and 24/tm maps do 
not overlap and we thus do not have 24/tm counterparts 



for matching. This decreases the effective area of the 70/im 
catalog from 11.5 to 10.5 sq. deg. 



4.2 Completeness Correction and Reliability Test 

We assess the completeness of the catalog as a function of 
flux density by injecting false sources at a given flux den- 
sity into the map and then determining what fraction are 
recovered by our source extraction techniques. An empirical 
noiseless point spread function image is used to add each of 
the point sources to the input map. The empirical PSFs were 
created using archival Spitzer data of 3C 273, 3C 279, and 
BL Lac. Data from individual AORs for these three galax- 
ies were processed into images using the MIPS DAT and 
then combined so as to filter out any extended, asymmetric 
emission. See Young et al. (2009) for additional details. 

We add 1000 artificial point sources in each of the flux 
density regimes where we test completeness, for a total of 
nine flux density bands and 9000 artificial sources. The re- 
sulting completeness as a function of source flux density is 
shown in Figure 7. Linear interpolation of this completeness 
function is then used to determine completeness correction 
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1.01620 


436+21+44 


2.64 


1.6 


1.3 


1.45 


650 


1.02180 


498+20+40 


2.70 


1.3 


1.0 


1.15 


1405 


1.02660 


606+16+61 


2.78 


1.0 


0.8 


0.90 


1907 


1.03100 


671+15+54 


2.82 


0.8 


0.6 


0.70 


4274 


1.03500 


806+12+8 


2.91 


0.6 


0.5 


0.55 


4179 


1.04700 


872+13+35 


2.94 


0.5 


0.40 


0.45 


7424 


1.06250 


952+11+48 


2.98 


0.4 


0.35 


0.375 


5785 


1.13600 


1006+13+30 


3.00 


0.35 


0.32 


0.335 


4424 


1.24000 


1051+16+84 


3.02 


0.32 


0.25 


0.285 


8140 


2.09150 


937+10+56 


2.97 



Table 1. ADF-S 24/im Number Counts. Columns give upper, lower and average flux density for each bin, the observed number of sources 
in this flux density range, the completeness correction factor for this flux density, the Euclidean normalized corrected counts and the log 
of these normalized counts. * indicates flux range where the cluster galaxies are likely to have increased the counts significantly above 
the field value. Two errors are given for the counts. The first of these is the Poisson error, the second is the LSS error estimated through 
a counts-in-cells analysis. The LSS error is given as where there are too few sources in the field to permit the counts-in-cells analysis. 
The two faintest data points correspond to the 80% and 50% completeness levels of the catalog respectively. 



factors in the calculation of the differential number counts. 
Wc find that our catalog is complete at the 80% and 50% 
levels at fluxes of 0.03 and 0.024 Jy respectively. 

We test reliability by using a negative map that other- 
wise matches the noise characteristics of the real map. This 
is produced in the same way as that used to test the reli- 
ability of the 24/im catalog. We find that our source list is 
close to 100% reliable over the full range of fluxes discussed 
here, with reliability falling from 100% to 99%, 97% and 
98% for the three faintest flux bins. This reliability assess- 
ment does not include the cross matching to 24/im sources 
so is a conservative estimate of the reliability of our cata- 
log. When we attempt to match the false, negative sources 
to genuine 24/im sources using the approach detailed above 
we find no matches at all. We thus conservatively conclude 
that our 70/im catalog is at least 97% reliable throughout 
the flux range. 



4.3 Photometric Correction 

A correction factor of 1.15 was derived by Frayer et al. (2009) 
in their analysis of COSMOS counts to account for the emis- 
sion outside the APEX PRFs. We use our injected artificial 
sources to test the correction factor needed in our data by 
examining the ratio of input flux density to recovered flux 
density. At bright flux density levels, comparable to the cal- 
ibration star HD180711 (S70 = 447.4mJy, Gordon et al., 
2007) used to check the Frayer et al. (2009) corrections, we 



measure the same correction factor of 1.15. At lower flux 
densities though, where the bulk of the 70/im counts lie, we 
find that a somewhat smaller correction factor is needed, 
as shown in Figure 8. We apply a linear interpolation of 
this correction function to our data in the calculation of the 
differential number counts. The overall absolute flux den- 
sity calibration of MIPS at 70/im is uncertain at the 5% 
level given that few of our sources are bright enough to be 
affected by nonlinearities in the detectors (Gordon et al., 
2007). 

4.4 70/im Counts 

The 70/im differential counts for the ADF-S are shown in 
Fig. 9 together with a variety of data and models from the 
literature. As with the 24/tm counts this is not a compre- 
hensive presentation of the full range of models currently 
available. The errors on the differential counts are a combi- 
nation of Poisson counting errors and LSS noise calculated 
using counts-in-cells. Our data set new constraints on the 
counts at bright flux density levels and place tighter con- 
straints at medium flux density levels than are currently 
available. They are in broad agreement with existing data 
where they overlap, though there is considerable scatter in 
all counts determinations at these wavelengths as a result of 
cosmic variance. At bright flux density levels we see a mod- 
erate excess of sources over what is expected on the basis of 
theoretical models. This may be a result of sources in the 
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Upper (Jy) Lower (Jy) Av. Flux density (Jy) Obs. Number Completeness Counts (gal sr 1 Jy 1 ' 5 ) log[(dN/dS)S 2 ' 5 ] 



2.00 


0.70 


1.35 


2.00000 


1.08050 


1104± 843 ±1766 


3.04* 


0.70 


0.40 


0.55 


9.00000 


1.04300 


2201± 765 ± 1321 


3.34* 


0.40 


0.23 


0.315 


22.0000 


1.03260 


2334 ± 514 ±933 


3.37* 


0.23 


0.14 


0.185 


29.0000 


1.01917 


1516 ± 287 ± 303 


3.18* 


0.14 


0.11 


0.125 


27.0000 


1.00917 


1573 ± 306 ±0 


3.20 


0.11 


0.09 


0.100 


36.0000 


1.00500 


1794± 300 ±0 


3.25 


0.09 


0.075 


0.0825 


39.0000 


1.00500 


1602 ± 258 ±0 


3.20 


0.075 


0.062 


0.0685 


61.0000 


1.00522 


1816 ± 234 ±0 


3.26 


0.062 


0.050 


0.056 


95.0000 


1.00710 


1855 ± 192 ± 


3.27 


0.050 


0.040 


0.045 


127.000 


1.01900 


1743 ± 158 ± 


3.24 


0.040 


0.030 


0.035 


190.000 


1.13000 


1543 ± 126 ±0 


3.19 


0.030 


0.024 


0.027 


165.000 


1.54800 


1599 ± 193 ±0 


3.20 



Table 2. ADF-S 70/im Number Counts. Columns give upper, lower and average flux density for each bin, the observed number of sources 
in this flux density range, the completeness correction factor for this flux density, the Euclidean normalized corrected counts and the log 
of these normalized counts. Errors on the counts are first the Poisson error and secondly the LSS error as determined by a counts-in-cclls 
analysis. Where the LSS error is the counts in cells analysis found no significant deviation in the counts in subficlds beyond what would 
be expected by Poisson statistics. * indicates flux range where the cluster galaxies are likely to have increased the counts significantly 
above the field value. The two faintest data points correspond to the 80% and 50% completeness levels of the catalog respectively. 
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Figure 7. 70/im Completeness function. Completeness as a func- 
tion of flux density for the 70/j.m observations. Completeness is 
determined by injecting artificial sources into the map and then 
determining what fraction are recovered by our source detection 
technique. A thousand artificial sources are used for the determi- 
nation of each point on the completeness function. 



foreground cluster at z=0.04 that lies in this field. This is 
also reflected in the large LSS errors found at bright flux 
levels. 



5 SOURCE CATALOGS 

The source catalogs for the ADF-S produced from this 
Spitzer data are available via the electronic version of this 
paper. The supplied data includes position and flux densi- 
ties for each source. Examples of the first lines of the catalog 
are given in tables 3 and 4 for the 24/xm and 70/xm catalogs 
respectively. 



Figure 8. 70/im Flux Density Correction function. Ratio of in- 
put flux density to output flux density for the artificial sources 
injected into the 70/im maps. This value is used to correct the 
photometry of the sources in our final catalog. 

6 DISCUSSION 

6.1 24-70pm Colours 

Each of our 70/xm sources is confirmed through cross iden- 
tification at 24pm, so we can immediately examine the 70- 
to-24 fim colour distribution for our sources. This is shown 
in Fig. 10. The colour distribution is consistent with that 
found by Frayer et al. (2006b) for the Spitzer FLS sources. 
If we assume a power law for the SEDs of our objects be- 
tween 24 and 70/im of the form f v oc v~ a we find our sources 
have an average spectral index of a— 2.8±0.3 compared to 
Frayer et al.'s value of 2.4±0.4. In the absence of redshifts, 
we cannot examine the detailed physics behind the varia- 
tion in colour. We note that sources with a log(70/24) <0.5 
are likely to have their far-IR emission powered purely by 
an AGN (Frayer et al., 2006b). We find 18 sources in our 
catalog with log(70/24)< 0.5 which can thus be regarded as 
candidate AGNs. However, one of these sources has such an 
exceptionally low ratio, log(70/24) <-l, that it is consistent 
with a Rayleigh- Jeans stellar continuum. We identify this 
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Figure 9. 70fim Euclidean Normalized Differential Counts. Euclidean Normalized Differential Counts for our 70/im data compared to 
other counts in the literature and to theoretical models. Data from this work is shown as solid dots. Literature counts are shown as 
follows: open squares: Frayer et al. (2006a) FLS warm+main field; open diamonds: Frayer et al. (2009) COSMOS field; open triangles: 
Frayer et al. (2006a) FLS verification field; stars: Frayer et al. (2006b) GOODS-N field; open circle: confusion analysis point from Frayer 
et al. (2006b). The models are plotted as lines, including: solid line: Rowan- Robinson (2009) model; dashed line: Pearson (2009) burst 
mode model; dash-dot line: Pearson (2009) bright mode evolution model; dash-three-dots line: Lagache et al (2004) model. The horizontal 
line above the counts indicates the range of flux densities expected for 0.2Ly () < L70 < 2L£ () objects in the z=0.04 cluster in this field 
and for which some enhancement in counts might be expected (see section 6.3). 



Number 


RA 


Dec 


Positional Error (deg] 





66.65229 


-55.20967 


0.0013 


1 


68.49929 


-55.04489 


0.0017 


2 


67.20997 


-54.50196 


0.0010 


3 


70.87384 


-52.54578 


0.0011 


4 


68.49582 


-54.16836 


0.0014 


5 


71.86896 


-52.33193 


0.0012 


6 


72.73022 


-53.46116 


0.0010 


7 


69.89899 


-53.01719 


0.0010 


8 


69.90396 


-53.01198 


0.0012 


9 


74.13552 


-53.02721 


0.0011 



Table 3. First ten lines of the 24/^m Spitzer ADF-S catalog 



24/im Flux Density (mjy) 24/im Flux Density Error (mjy) 



284.95 


0.067 


271.58 


0.020 


188.45 


0.020 


173.64 


0.020 


141.12 


0.020 


115.68 


0.059 


115.21 


0.060 


110.85 


0.071 


109.52 


0.056 


101.69 


0.020 



source as a Doradus, a bright binary star system (V=3.25 
and V=4.3) whose primary is a chemically peculiar vari- 
able star (Heck et al., 1987). It should be noted that this 
source is successfully identified as a star through our K-[24] 
colour analysis and is thus not included in our considera- 
tion of the galaxy counts. We thus conclude that 17 of our 
sources have their mid-to-far-IR luminosity powered by an 



AGN. This fraction of AGN sources is consistent with the 
FLS results. Figure 10 also shows the range of colours that 
would be found for three archetypical objects over a range 
of redshifts where they might be detected. These include the 
local, quiescent spiral galaxy M100, shown from z=0 to z=l, 
a generic AGN SED, the local starburst M82, and the near- 
est ULIRG, Arp220, all shown from z=0-5. This suggests 
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M82 

AGN""'" Arp220 
M100 



Log(F70/F24) 

Figure 10. Histogram of the 70 to 24 fim flux density ratio for 
the sources in the 70/^m catalog. The colour distribution is sim- 
ilar to that seen by Frayer et al. (2006b) in the FLS. The single 
source with a very low 70/24 ratio is a bright binary star system, 
a Doradus whose flux density ratio is consistent with its 70 and 
24 (im emission coming solely from the stellar photospheres. Also 
shown are the ranges of colour that would be found for various 
sources over a redshift range where they are likely to be detected. 
These include a generic AGN SED, the starburst M82 and the 
nearest ULIRG, Arp220, shown for z=0-5, and the lower lumi- 
nosity quiescent local spiral galaxy M100 shown for z=0-l. Note 
that the colour ranges are not monotonic with redshift so, for 
example, the reddest colour for Arp220 occurs at z~1.5. 



that many of the sources detected at both 70 and 24/xm in 
this sample are actively star forming objects like M82 or 
Arp220. 



6.2 Count Models 

We compare our observed number counts, and those from 
the literature, to a number of different count models. These 
include models from Lagache et al. (2004) and Rowan- 
Robinson (2009) as well as two models from Pearson (Pear- 
son & Knan, 2009) which have yet to be fully described in 
the literature. 



6.2.1 Pearson Count Models 

The Pearson model incorporates an infrared backward evo- 
lution framework following the models of Pearson (2005) 
and Pearson et al. (2007). These models were previously 
successfully used to reproduce the combined mid-infrared 
source counts from ISO & Spitzer at 15pim & 24/im and 
have recently been successfully applied to the source counts 
at submillimetre wavelengths (Pearson & Khan, 2009) The 
model uses the 60/xm luminosity function derived from the 
IRAS Point Source Catalogue (Saunders et al., 2000) to rep- 
resent cool and warm galaxy populations defined by local 
IRAS colours, where cool 100/im/60/im cirrus-like colours 
represent the quiescent, non-evolving normal galaxy popu- 
lation and warmer 100/im/60/im colours represent evolving 
star-forming galaxies. To model the AGN, the luminosity 
function derived from the selected sample of Matute 

et al. (2006) is used. The model incorporates a suite of spec- 
tral energy distributions to model the cool quiescent, warm 



star-forming and AGN populations. The spectral templates 
for the quiescent normal galaxy population have been se- 
lected from the libraries described in Efstathiou & Rowan- 
Robinson (2003). The adopted star- forming templates are 
selected as a function of increasing luminosity for starburst 
(L IR < 10 n L Q ), LIRG (lO n L < L m < 1O 12 L ) and 
ULIRG (L IR > 1O 12 L ) populations. A total of 7 SEDs (2 
starburst, 3 LIRG, 2 ULIRG) are selected from the starburst 
model template libraries of Efstathiou et al. (2000) to ensure 
a reasonable variation in the mid-infrared PAH features in 
an attempt to avoid features in the source counts caused by 
a single choice of template for all sources. The AGN tem- 
plate is taken from the tapered disc dust torus models of 
Efstathiou & Rowan-Robinson (1995). All the spectral tem- 
plates have been shown to provide good fits to IRAS, ISO 
and Spitzer selected samples of infrared populations (Efs- 
tathiou et al. (2000) , Rowan- Robinson et al.(2004), Rowan- 
Robinson et al. (2005)). 

The general Pearson Model framework is available in 
two evolutionary flavours referred to as the bright and burst 
models categorized by their dominant populations of star- 
burst (M82 like) and LIRG/ULIG (Arp 220 like) sources 
respectively. The bright model broadly follows the evolution- 
ary scenario of Pearson & Rowan-Robinson (1996) and as- 
sumes all active populations evolve in both luminosity and 
density as a double power law, of evolutionary strength k of 
the form (l + z) k . The evolution rises steeply to a redshift of 
unity and more shallowly thereafter. The burst model was 
originally introduced in Pearson (2001) and assumes that the 
more luminous LIRG & ULIRG populations evolve exponen- 



tially (of the form 1 + k.exp[ 



], where k and a are 



evolutionary parameters) from the present epoch to some 
peak redshift z p ensuring both local rarity and the emer- 
gence of the population towards a peak redshift ~1. The 
bright and burst models are explained in detail in Pearson 
(2009). 



6.2.2 Comparison to Observed Counts 

At 24/im all of the models perform reasonably well with no 
gross disagreements with the observations. However, none 
of the models are an ideal fit to the data. The Pearson 
burst model underpredicts the counts below ~0.3mJy while 
the other three models do well in this flux density regime. 
Between 0.3 and ~3mJy all models except the Rowan- 
Robinson model match the data well, with the Rowan- 
Robinson model overpredicting the counts around l-2mjy. 
Counts at flux densities brighter than a few mjy are still 
quite poorly constrained, especially considering the problem 
with the cluster in the current data set. The best constraints 
come from SWIRE and suggest that all models except that 
of Rowan-Robinson are overpredicting the brightest counts. 
The Lagache model has the most significant disagreement 
in this respect. The all sky WISE survey (Eisenhardt et al., 
2009) will settle the issue of the brighter counts now that 
this satellite has been launched. This visual inspection of the 
counts marginally favours the Pearson bright model, with 
the Lagache model second followed by the Rowan- Robinson 
and Pearson burst model. 

At 70/im, possibly because this is a wavelength regime 
that was not well studied prior to Spitzer, there are rather 
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Source Number RA Dec Positional Error (deg) 70/im Flux Density (Jy) Signal/Noise ratio 



0.0 66.96191 -55.05889 0.00068 0.92 125.0 

1.0 68.49582 -54.16836 0.0014 0.91 134.6 

2.0 68.81533 -54.31517 0.0009 0.80 122.7 

3.0 69.10008 -54.41525 0.0008 0.68 111.2 

4.0 68.43724 -53.77945 0.00087 0.66 95.2 

5.0 69.54994 -54.37454 0.00054 0.60 94.4 

6.0 68.59667 -54.69228 0.0013 0.60 96.6 

7.0 70.02633 -54.34266 0.00058 0.59 86.1 

8.0 73.56086 -53.36012 0.00069 0.59 83.7 

9.0 66.76700 -54.26271 0.00059 0.56 78.9 



Table 4. First ten lines of the 70/im Spitzer ADF-S catalog. The positions given in the 70/j.m catalog are those of the associated 24/im 
sources since the 24/im positions are more accurate than those obtained at 70pm. 



larger differences among the models and between the data 
and the models. The Rowan-Robinson model significantly 
underpredicts counts in the 0.02 to 1 Jy range and fails to 
place the peak of the normalized counts at the right flux 
denstiy. The Pearson burst model overpredicts counts in the 
0.002 to 0.01 Jy range, predicting too broad a peak for the 
counts. The Pearson bright and Lagache models both visu- 
ally match the data reasonably well. 

We perform a quantitative test of the models against the 
data by calculating the sum of the x 2 values of the models 
compared to the data sets shown in Figures 6 and 9. These 
support the conclusions drawn from the visual inspections of 
the data. At 24pm the Pearson bright model has the lowest 
X 2 by a considerable margin, with less than half the % 2 value 
of the Rowan-Robinson model and less than a third the \ 2 
of the Lagache or Pearson burst models. At 70pm the re- 
sults are less clear cut, with the Lagache and Pearson burst 
and bright models having very similar \ 2 values while the 
Rowan-Robinson model has roughly three times this value. 

We conclude that the Pearson bright model fits the 
counts at 24 and 70pm more accurately than the other mod- 
els examined here. Additional data at the bright end of the 
24 pm counts, at the fainter end of the 70pm counts, and flux 
densities close to the peak of the 70pm counts (i.e. at 0.01 
Jy and fainter) would be the most useful for constraining 
the models further. There is currently only one Spitzer data 
set providing constraints to the 70pm counts at O.OUy and 
fainter, the GOODS-N counts from Frayer et al. (2006b), so 
more data in this regime are particularly needed. Such ob- 
servations should be forthcoming from both the Spitzer Far- 
Infrared Deep Extragalactic Legacy (FIDEL) survey (Dick- 
inson et al., 2007) and from observations with the PACS 
instrument on the Herschel Space Observatory. The latter 
is also able to probe deeper in small fields because of the 
fainter confusion limit provided by the 3.5m Herschel pri- 
mary mirror, as demonstrated by initial results at 100 and 
160pm from the PEP survey (Berta et al., 2010). 



6.3 The Local Cluster - DC 0428-43 

As discussed in the introduction, the ADF-S encompasses 
the z=0.04 galaxy cluster DC 0428-53 (Dressier, 1980; also 
known as Abell S0463). We might thus expect some en- 
hancement in the counts seen at flux densities appropriate 
to ~ L* galaxies at the cluster redshift. We take estimates 
for the field galaxy population L* at 24pm from Babbedge 



et al. (2006) and for 70pm from Patel et al. (in preparation), 
yielding flux densities of 5.8mJy and 310mJy respectively, 
and show an indicative range of effect from 0.2 — 2 L* in 
Figures 3 and 6. As can be seen there is some indication 
of enhanced counts at the brighter end of this flux density 
range at both 24 and 70pm with less of an effect at fainter 
flux densities. This is to be expected as the contribution from 
cluster objects will be less significant at fainter flux densities 
where we detect greater numbers of objects in the field pop- 
ulation. For example, in the 3.5mJy flux density bin at 24pm 
we detect 290 objects while Dressier et al. (1980) detects a 
total of only 131 optically selected galaxies in the cluster 
with 84 spectroscopically confirmed (Dressier & Schectman, 
1988). There is some indication of excess counts above model 
expectations or SWIRE observations in the brightest 24pm 
bins, equivalent to 10 — 20 L| 4 at the cluster redshift. Cluster 
objects might be responsible for this if their 24pm luminos- 
ity is enhanced beyond that of the field population, possibly 
through enhanced star formation rates resulting from in- 
teractions between galaxies in different cluster components, 
and we do indeed find that four of our five brightest objects 
lie in the third of the field adjacent to the cluster. Such ef- 
fects have been suggested elsewhere (Fadda et al., 2000). In 
this context it is interesting to note that Dressier & Schect- 
man (1988) find evidence for substructure in DC 0428-53, 
with a colder subsystem apparently merging with the main 
cluster. 

We attempt to remove the influence of the cluster from 
our counts by excluding from consideration those sources in 
the 24 and 70pm catalogs that are coincident with sources 
in the Dressier & Schectman (1988) catalog of objects spec- 
troscopically identified as cluster members. This is shown 
in Figs. 6 and 9. As can be seen this removes the enhance- 
ment at bright fluxes seen at 70pm, but does not have the 
same effect at 24pm. This would suggest that there are clus- 
ter members detected at 24pm that are not in the list of 
spectroscopically confirmed cluster members. Given that the 
24pm observations are quite deep and that the optical spec- 
troscopy is quite shallow, V<16, this is not unexpected. Fur- 
ther examination of this cluster and its galaxies is underway 
but is beyond the scope of this paper. 



7 CONCLUSIONS 

We present the results of 24 and 70pm Spitzer observations 
of the ADF-S region, which covers a contiguous field of ~12 
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sq. deg.. We derive galaxy number counts from this data. 
They confirm that the galaxy population in these bands 
evolves strongly. We compare our own and literature counts 
to galaxy count models from Lagache et al. (2004), Rowan- 
Robinson et al. (2009) and two models by Pearson (Pearson 
& Khan, 2009). We find that the Pearson bright model pro- J 
vides a reasonable fit to the counts in both bands, perform- o 
ing as well as, or better than, the more established Lagache ~ 
and Rowan-Robinson models, but the Pearson burst model S, 
does not perform as well. None of the three best perform- <g 
ing models provide a fully successful match to the data. We £ 
also discuss the contribution of galaxies within a local clus- " 
ter, DC0428-43 at z=0.04 which lies within the ADF-S, to 
the counts and find indications that some cluster galaxies 
may have their 24/im luminosities enhanced beyond what is 
found in the field population. The full source catalogs de- 
rived from these observations are available in electronic form 
for download. 
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APPENDIX A: STELLAR COUNTS AT 24 
MICRONS 

In Figure Al we provide a plot of the number counts of 
sources identified as stars by our K-band vs. 24/nm colour 
analysis. 



